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Abstract 

In this paper we study a storage process or a liquid queue in which 

^ ' the input process is the local time of a positively recurrent stationary 

OO . diffusion in stationary state and the potential output takes place with 

zL ! a constant deterministic rate. For this storage process we find its sta- 

\^ I tionary distribution and compute the joint distribution of the starting 

Cp . and ending times of the busy and idle periods. This work completes 

^—N , and extends to a more general setting the results in Mannersalo, Nor- 

■^ \ ros, and Salminen j^]. 
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1 Introduction 

Let A = {At : t G M} and B = {Bt : t G M} be increasing stochastic 
processes. We view A as the input process to a storage or a buffer and B as 
the potential output. One of the standard ways to define a storage process 
associated to A and B is via Reich's formula (see Reich [TT] ) 

St:= sup {At-A,-{Bt-B,)}, teR. 

—oo<s<t 

We refer also to Prabhu jTU] P- 114 for an approach in which the storage 
process is defined as a solution of a stochastic integral equation and the 
input process A is a subordinator without drift and Bf = t. In Harrison [7] 
p. 14 the case in which A is a Brownian motion and Bt = fit is studied via 
Skorohod's reflection equation. Brownian motion, although not increasing, 
is much analyzed input process and its appearence in these models can be 
motivated in different ways, see, e.g., Harrison [7j p. 30, Roberts, Mocci and 
Virtamo ^2] P- 377, and Salminen and Norros [?]. A new kind of a model is 
introduced in Mannersalo, Norros, and Salminen pi, in which the local time 
at of a reflected Brownian motion with negative drift serves as the input 
process. Recall that this local time process is increasing and continuous but 
the points of its increase form a singular set. 

In this paper, we consider a model of a storage or a fluid queue, where the 
input is the local time L at zero of an arbitrary one-dimensional positively 
recurrent stationary diffusion and the potential output Bt = fit, fi > 0; thus 
extending the approach in 0|. Recall that any subordinator can be viewed 
as the inverse local time of a strong Markov process. Hence, in a sense, the 
local time storage models are companions of the subordinator storage models 
and as such interesting objects of research. 

It is shown in [9^ that the storage model studied therein can be obtained 
via a limiting procedure from some on/off processes. Especially, a priority 
discrete queueing system is constructed which, when in heavy traffic, leads 
to a storage process with local time input. The proof of weak convergence 
in fU' is based on the continuous mapping theorem and the interpretation 
of the Brownian local time as a supremum process. However, this approach 
is not applicable in our general case, and we do not discuss here the weak 
converegence aspects of our model. 

Further, in [H] the joint distribution of the starting time gi and the ending 
time di of an idle period and the marginal distributions of the starting time gi, 



and the ending time df, of a busy period are obtained. The joint distribution 
of gb and d^ is not, however, given in p. The aim of this paper is to fill this 
gap and to generalize the results to an arbitrary diffusion local time storage. 

We give explicit expressions for the joint Laplace transforms of the start- 
ing and the ending times of the busy and idle periods. Surprisingly, the 
answers have a short simple form, expressed in terms of the Levy-Khinchin 
exponent of the inverse of the local time process L. By considering a marked 
point process obtained from the storage process S we make the simple but 
crucial observation needed for finding the joint distributions. This observa- 
tion, formulated here only for the busy period, is that the random variable 
(— (?h, d}f) is identical in law with {UV^ (1 — U)V)^ where U has the uniform 
distribution on (0, 1) and V^ is a positive random variable independent of U 
having the distribution of the length of the busy period. In particular, which 
is important here, we show that the joint Laplace transform of such variables 
can be easily expressed in terms of the marginal Laplace transforms. 

We remark also that in the marginal distribution of —gb is computed 
by making use of Bertoin's path decomposition theorem (see Bertoin j2j). 
However, in this paper we focus instead on di,- It is seen that the formula for 
the first time when a spectrally positive Levy process jumps over a level, see 
()3.9p . can be applied to find the Laplace transform of dh- 

The paper is organized as follows. In the next section we give some 
definitions and construct the storage process with a diffusion local time input. 
It is important to verify here that our local time process has stationary 
increments. It is also proved that the stationary distribution of this storage 
process is always exponential with an atom at zero - only the value of the 
parameter and the size of the atom vary when the underlying diffusion is 
changed. In Sectional the Laplace transform of the ending time di, of a busy 
period observed at time zero is computed. In Section |5 we find the Laplace 
transform of the starting time gi of an idle period observed at time zero. 
Theorem 15.11 in Section El can be viewed as the main result of the paper, 
where the joint distributions of {—gb,db) and {—gi,di) are given via their 
Laplace transforms. Finally, in Section [HI we discuss the case treated in p. 
In particular, we invert the Laplace transforms to obtain the joint densities 
of {-gi,di) and {-gb,db). 



2 Definitions and preliminaries 

Let X = {Xt : t > 0} be a one-dimensional time-homogeneous, regular, 
conservative diffusion living in an interval / C R. We assume, without loss 
of generality, that G /. The probability measure and the expectation 
associated with X started at x, are denoted by P^ and E^^, respectively. 
Recall that X has a symmetric transition density with respect to its speed 
measure m: 

PxiXt e dy) = p{t; X, y)m{dy) = p{t; y, x)m{dy). 

We introduce also the Green function 

/■oo 

Ga{x,y)= e''^^p{t;x,y)dt, x,yel,a>0. 
Jo 

Further, it is assumed that X is positively recurrent, that is, the speed mea- 
sure is finite in I. Denote M := m{I} < +cxd. For more details on linear 
diffusions we refer to Ito and McKean jS] and Borodin and Salminen 0. 

Let X(i) and X^^) be two copies of X. We take X^^^ = xf'^ with the 
common distribution 

P(xi') G dx) = P(Xi') G dx) = ^^^ =: m{dx) 

but otherwise we let X*^^^ and X*^^) be independent. We remark that fh is 
the stationary probability distribution of X. Define Z = {Zt : t G M} as 



Zt-.-- 



X3 if t < 0, 
Xf ^ if t > 0. 



The process Z thus defined is a stationary process in stationary state. 

Let {Ll : t > 0} and {Ll : t > 0} be the local times at zero of the 
processes X*-^-* and X*-^-*, respectively. We choose the normalization of L*-'-*, 
i = 1,2, with respect to the speed measure, that is, 

L« := hm ./ .. / l(-e,.)(X«) ds, z = 1, 2. (2.1) 

Using L^^) and L^'^\ we define the local time process L = {Lt : t G M} via 



^(1) 

-(2)* 



if t < 0, 



L^^ if t > 0. 



The process L is an increasing continuous process passing through zero at 
time zero. Next we give an important property of L, which allows us to 
construct a stationary storage process. 

Proposition 2.1. The process L has stationary increments, i.e., for all 
t > s, 

Lt — Ls = Lt-s- 

Moreover, 

E(Li) = t/M. (2.2) 

Proof: Consider first the case i > s > 0. Then, using that Z is a stationary 
process in stationary state, we have 



P{Lt-Lsedl) = m{dx)P{Lt-L,edl\Z, = x) 

= / fh{dx)Px{Lt-s e dl) 
= P(Li_, Grf/). 

The case s < t < is treated similarly. Assume next that s < < t. Then 
by the conditional independence of X^^^ and X^"^^ given Zq we can write 



-aL^'h 



= /"m(rfx)E,(e-"^*'')E,(e- 
On the other hand, consider 

= fm{dx) /"E,(e-"^- ; Z_, G dy)E(e"(^*--^-) | Z^, = y). 

Let now r be an exponentially distributed random variable with parameter 
a independent of X^'^\ and define the killed process 

;,.._/xf if Lr<r, 
^•"l A ifLf)>r, 

where A is a cemetery point. Then (see Borodin and Salminen [S] p. 28 
or Ito and McKean p. 179-183) X* is a diffusion having the same speed 



measure as X^^\ Let p*{t;x,y) be its (symmetric) transition density with 
respect to m. We clearly have 

P^iX^edy) = P.(xf edy,Lf <r) 

= P*it]x,y)m{dy). 
From the symmetry of p*(t; x, y) it follows that 

/"m(rfx)E,(e-°^- ; Z_, G rfy) = m(t/y)E,(e-"^-0, 

and, therefore, 

= /"m(rfy)E,(e-"^-OE,(e-"^') 
= E(e""(^*"^=)), 

as claimed. Finally, to show (j2.2|) recall the formula (see Ito and McKean jH] 
p. 175 or Borodin and Salminen [^ p. 21) 

E^(Lt) = / p(s;x,0)ds, t > 0. 
Jo 

Consequently, 



because 



E(Li) = m{dx) / (isp(s;a;,0) 

(is / ■m{dx)p{s; x, 0) 
'0 i/ 

ds / m{dx)p{s] 0,x) 
/o A 

t/M, 



■m{dx) p{s; 0,x) = 1. 



D 



Definition 2.2. The process 5* = {S't : t G M} defined as 

St := sup {Lt — Ls — fi{t — s)}, t eM. 

—ao<s<t 

is called a storage process with local time input, constant service rate fj,, and 
unbounded buffer associated with the process Z . 

Because the process L has stationary increments it follows from Definition 
12.21 that S" is a stationary process (in stationary state). Our next task is to 
find its stationary distribution, that is, the distribution of 5*0. Notice that 
5*0 is determined by the process X^^\ With a slight and short abuse of 
our notations, let {Lj : t > 0} be the local time of X^^^ started at zero 
normalized as in ^Z.\\ . Then the process A'^^' given by 

Af ^ := inf{s : L^ > t}, t > 

is a subordinator (a right-continuous increasing process on M_|_ with indepen- 
dent and stationary increments). When normalizing as in (|2.H) we have (see 
Ito and McKean 0, p. 214) that 

E„(exp{-„^<')})=exp|-^}, (2.3) 

and, therefore, 

Eo (exp{-«(Af ) - t//x)}) = exp l-t ( ^ ". " ^) } = e^P{^^(«)}' 
where the function 

*(") - ^ Glib) ''■^' 

is called the Levy-Khintchin exponent of the spectrally positive Levy process 
{A«-t//i : t>0}. 

Proposition 2.3. Assume that fi > 1/M, in other words fi > E(Li). Then 
the equation 

^p{a) = (2.5) 

has a unique positive solution a*. Moreover, 

1) P{So > t\Zo = x) = E (e-°*^o|Zo = x) e""^*, t > 0, x G /, 



where Hq := inf{t > : Zt = 0}; 

2) the process S is stationary with the stationary distribution given by 

P(^* > y) = ^e-V^ y>0, 

andP{St = 0) = l~j^. 

Proof: The stationarity follows from Proposition l2.1[ The other claims follow 
from the fact that the spectrally positive Levy process {A\ — t/fi : t > 
0} drifts to +00 and its infimum is exponentially distributed in M_ with 
parameter a*. For further details, see [Tn|. □ 

Remark 2.4. 1) The Laplace transform of Hq appearing above can be ex- 
pressed as 

*^^^' ^"G«(0,0)' ^ ^^^ 

cf. Borodin and Salminen j3], p. 18. 

2) Recall (cf. Bingham jl] p. 720) that for a > a*, the function ip{a) 
is positive and increasing and, hence, has an inverse which we denote by 
a f— i> ri{a). 

3) From Proposition 12.31 it is clear that So, conditionally on So > 0, is 
exponentially distributed with parameter a*/fi. 

We assume throughout the paper that ^ > 1/M. Next we define the 
subject of our main interest in this work, the on-going busy and idle periods 
of ^. 

Definition 2.5. If 5*0 = then the random variables 

Qi = sup{t < : St > 0} and di = inf{t > : ^t > 0} 

are called the starting time and the ending time, respectively, of the on-going 
idle period at time zero. If 5*0 > then the random variables 

gb = sup{t < : Si = 0} and di, = inf {t > : S^ = 0} 

are called the starting time and the ending time, respectively, of the on-going 
busy period at time zero. 



3 Busy periods 

In this section we compute the Laplace transform of the ending time of the 
on-going busy period. It is shown in Mannersalo, Norros, and Salminen jH] 
using a point process view on the busy and idle periods that —gb and df, have 
the same law. In the present paper (in Section Ej) we develop this approach 
in order to find the joint distribution of —Qb and db- 

From Proposition 12.31 it is seen that the probability that there is a busy 
period at time zero is l/(M/i). 

Theorem 3.1. Lett] be the inverse of tp (cf. Re'mark \2.4\ 2). Then 

E(e-"'^^ I So>0) = -^y (3.1) 

The rest of this section is devoted to the proof of Theorem 13. II Assuming 
that there is a busy period at we have (cf. pi) 

St = Lt - fit + So, gb <t < db. 

Let Hq be the first hitting time of zero for the process X^'^\ Because 



(1)1 

s<0 

r{2) 



5*0 = sup{/is — Lg} = sup{/is + L_l} 

s<0 s<0 



it follows that 5*0 and Hq are conditionally independent, given Zq . The 

(2) 

local time Lj = when < t < H^ . Hence, 

St = So — fit, < t < minjifg ; So/fi}- 
To proceed consider, therefore, two cases: 

1) So/ft < H^o\ 

2) So/fi > H^o^ (see Figured). 
In the case 1) we obviously have 

[So/fi < H^o^ , 5o > 0} = {db = So/fi , So > 0}. (3.2) 

To analyze the case 2), we study in detail some relationships between So and 

(2) 

Hq '. Firstly, introduce 

J{a,p) := /"m(t/x)E^(e-°-^°)E^(e-^^°), a > 0, /? > 0. (3.3) 

8 



Proposition 3.2. 



E (e-(^o/.-H<^V/3Hf > , s,/^, > Hl^'A = ^-^J{a\ «* + /?). (3.4) 
V / a + a* M 

Proof: Using Proposition 12.31 and ttie conditional independence of 5*0 and 
Hq given Zq, we compute the joint Laplace transform as follows: 

//•OO 
midx) / E,. (e-"(^o^-t)e-/3t ^ ^^/^ > ^^ p^(^(2) ^ ^^) 

/• /"OO /"CO 

= / m{dx) / e("-^)*P^(i7(5^^ G rft) / e-°"P^(So//i G du) 
J I Jo Jt 

= / m(dx) / e("-^)*P,(i7i'^ G dt) / e-""a*e-"*"E,(e-"*^°)ciM 

— fm{dx) f e("-^)*P,(i7i') G rft)E,(e-°*^°)e-("+"*)* 
M J I Jo 



a + a* M 
a' 



1: y"m(rfx)E.(e-°*^«)E.(e-(-*+^)^«; 



a + a* M 

as claimed. D 

An immediate consequence of Proposition 13.21 is the following 

Corollary 3.3. 1) The random variables S'o//i — Hq and Hq are condi- 
tionally independent, given that So/ ft > Hq . 

2) The conditional distribution of Sq/^ — Hq , given that So/n > Hq , is 
exponential with parameter a*. 

f,?> o /... rr{2)\ J{a\a* + a) 



3) E (e~"^o^' , So/fi > Hf 

4) E ('e-"^o/'^ , 5'o//i > //, 



M 
a* J{a*,a* + a) 



(2) 

° / a + a* M 

Proposition 3.4. The functions J {a*, a* + a) andipla) defined in ()c/.5|) and 
( li^.y^D , respectively, satisfy the relationship: 

T/ * * N 1 ip{a + a*) .^ ^. 

J{a*,a* + a) = — -. (3.5) 

H a 



Proof: Using formula (j2.6|) and the Chapman- Kolmogorov equation, we 
compute for any positive ai 7^ 02, 

m{dx) E^(e-"^^°)E^(e-"2^«) 

midx) G^, (x, Q)G^, (x, 0) /{G^, (0, 0)G«, (0, 0)) 
I 

t*co 1*00 

m{dx) / dte-'''^p{t;x,0) / dse-'^''p{s;x,0)/{Ga,{0,0)Ga2{0,0)) 
I Jo Jo 

00 /"OO /• 

rft / dse-"^*-"2" / m(rfa;)p(t;x,0)p(s;x,0)/(G'«,(0, 0)^^2(0,0)) 
Jo Ji 

/•oo /»oo 

dt rfse-"i*-"^"p(t + s,0,0)/(Ga,(0,0)G«,(0,0)) 
Jo io 

rfw ( f e("^-°^)'rfsJe-"^XM;0,0)/(Ga,(0,0)Ga2(0,0)) 

1 /""^ 

'g-«2« _ e-"i")p(M; 0, 0)rfM /(G^i(0, 0)Ga,iO, 0)) 



"1 - "2 Jo 

1/1 1 



«i-a2 VG'«i(0,0) G'a2(0,0) 
Since 



a 



G.(0,0) 




u:j, 


J(ai,a2) = - - 
/i 


ai - 


-^(aa) 
- 02 



we have 

1 ih(n/-,\ — ihln/r,) 

(3.6) 

Substituting a* for ai and a + a* for 0^2 in ()3.6p and using that ip{a*) = 
give (jSS)). D 

To proceed with the proof of Theorem 13. II write 

E(e~"'^'' I 5o > 0) = E{e-'"^' , So/fi < H^^ | 5o > 0) 

+ E(e-"'^SSo//">i7i')|5o>0). (3.7) 

For the first term we obtain by ()3.2p . Corollary 13.31 4. and Proposition 12.31 

E(e-"^'' , ^o//i < H^o^ I 5o > 0) 

= E(e-"^°/'^ , ^0^ < H^'^ I So > 0) 

10 



E(g-aSo^ I 5g > 0) - E(e-"^o/^ , ^o//i > H^^^ \ So > 0) 



E(g-«So^ I 5-^ > 0) - 



P{So > 0) 



a' 



a + a 



-(l — fiJ{a*, a* + a)). 



(3.8) 



To compute the second term in (j3.7p consider the process 

^ = {^H(^)+. ^ ^ > 0}- 

By the strong Markov property X has the same law as X^^' started at zero. 
Let L = {Lt : t > 0} be the local time at zero of X, that is, the local time of 
X*^^^ started at zero. Let A be the right continuous inverse of L and introduce 

T+ := inf{s : A, - s/fi >t}, t> 0. 
Observe next that in the case Sq/^ > Hq (see Figure [T} we have 

"fe - -^c. , „(2) H • 

Using this representation and conditional independence stated in Corollary 
13.31 1. it is seen that 



E (e-'^^ , So/fi > Hf^ I So > 

= E (e-"'^' I 5o//i > i^r) P [So/ii > Hf^ I 5o > 
= Efexp|-af-T+ ,,, + ^ 
xP (5o//i > H^'^ I 5o > 



/i 



5o/^ > H'o 



(2) 



E exp <^ -a - T 



r(2) 






(2) 



5*0 //i > H^ 



(2) 






-a I — T^ „^ H Hr, 



(2) 



/^ 



E I exp 

xE (e-°^o'' , S'o/^ > H^^) / P (So > 0) . 



5o//i > if, 



(2) 
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From Corollary 13 .31 2 and the fact that T"*" is independent of 5*0 and Hq , we 
obtain 

= /" E (exp I --r+ - atl j P (So/fx - H^^^ e dt \ So/fi > H^^'^) 

= a* He (exp \--tA) e-("+"*)*dt. 

To compute the last integral we recall the following formula for the double 
Laplace transform of the first time when the spectrally positive Levy process 
{As — s/fi : s > 0} jumps over the level t > (see Bingham jl|, p. 732): 

where a ^-^ '?/'(«) is given by p.4|l and a ^-^ ^(«) is the inverse oiip{a), a > a*. 
Combining this with the results of Proposition 12.31 and Corollary 13.31 yields 

E (e-"'^" , So/fx > H^o^ I ^0 > O) 

a* I f (^(f) -(« + «*)) \ 

= -— — 1 - ^ ^ / V ^^ Jioi*^ «* + «). (3.10) 

" + " V ^(f)(f-^(« + «*))/ 

Finally, recalling ()3.7|) . putting ()3.10|) together with ()3.8p . and using ()3.5p 
complete the proof of Theorem 13.11 
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Figure 1: The on-going busy period at time zero. 
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4 Idle periods 

Let Hq be the first hitting time of zero for the process X^^'^ and introduce 

^ := sup{/it - L_ m }. 

By the strong Markov property of X^^^ and Proposition 12.31 1. the random 
variable ^ is exponentially distributed with parameter a*/n (cf. Remark 
12.41 3) and independent of ifg . The local time Lt is equal to zero when 

^0 



-ifn < t < 0, and for these values of t 



St = supjyus - Ls+t}- 

s<0 

Since ^ = >S'_„{i) > it is clear that —Hq < gi < 0, a.s displayed in Figure 
121 Further, we have (see Figure |2I) that 

{So = 0} = {e/M < 4'^}, 

and that the conditional distribution of —Qi, given that Sq = 0, is the same 
as the conditional distribution oi Hq — ^//x, given that ^//i < Hq . 

Theorem 4.1. For a ^ a* , 

a [Mjj, — l)(a — a*) 
Proof: Using the independence of ^//i and Hq and ()3.6p . we compute 
E(e"^' , So = 0) 

p poo 

= / m{dx) / P.(iy,S'^ e rft)e-"*Eo(e°«/'^ , ^/^i < t) 
Ji Jo 

= I m{dx) I V^iHi^^ E dt)e 



oo ,, 



a* — a 



a 



M{a* - a) J J 

14 



y"m(rfx)(E.(e-"^°)-E.(e-"*^«)) 



a* fl i){a) 1 i){a*) 



M{a* — a) \fi a fi a* 



Ma{a-a*)' 
Noting that 



P(^o = 0) = 1 ^ 



gives dHH). n 

Remark 4.2. In case a = a*, (|4.ip should be read as 



15 




Figure 2: The on-going idle period at time zero. 
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5 The joint distribution of the starting and 
the ending times of the on-going busy and 
idle periods 

In this section, as our main result, we compute the joint distribution of the 
starting time and the ending time of the busy and idle periods. We remark 
that it is possible to compute the joint distribution directly in the case of the 
idle period, but we have not been able to do this in the busy period case. To 
resolve this difficulty we use the theory of Palm probability and properties of 
a special class of two-dimensional random variables to show that it is enough 
to know only the Laplace transforms of db (respectively, —gi) in order to 
obtain the joint Laplace transforms of db and —gb (respectively, di and —gi). 

Theorem 5.1. 1) For a ^ P, 

K(e«.-|.o>0)^^(^-^). (5.1, 

2) For ay^ (3, 

^(^ "^° = °) = (M^-l)(«-/3)l^3^-^3^J- ('•') 

Remark 5.2. Letting /5 — > a in (I5.1|l and ()5.2|1 gives us the Laplace trans- 
forms of the lengths of the busy and idle periods 

E (e-C-^.-..) I ^„ > 0) = «*-^ f-^) 
^ da \ri[a/jj,)J 



and 



respectively. 



^ ' ^ M^-lda \a-a*J 



To prove Theorem 15.11 we use the point process approach to the storage 
process 5" proposed in Mannersalo, Norros, and Salminen |9^. When S hits 
zero it stays there for a positive amount of time a.s. Hence, we can construct 
from 5* a stationary marked point process N = {{Tn, Zn) : n G Z} with a 
finite intensity A, where for every n, Tn and T„+i are the starting and ending 
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times, respectively, of an idle or a busy period and Z„ is the mark associated 
with point T„. We let Z„ = if T„ is a starting point for an idle period, 
otherwise Z„ = 1. Let Pn and JF denote the probability measure and the 
natural a-algebra, respectively, induced by A^. We assume that {T„ : n E Z} 
are ordered so that — oo < ■ ■ ■ < T_i < Tq < < Ti < T2 < ■ ■ ■ + oo. Let 9t 
be the shift operator and P^ the Palm probability, which is related to the 
probability measure P^ via Ryll-Nardzewski and Slivnyak's formula: 

POO 

Piv(A) = A / P%{Ti >t,ete A) dt, AeT. (5.3) 

Recall that the Palm probability Pat can be interpreted as the conditional 
probability, given that there is a point at time zero, i.e., 

P^d^o = 0}) = L (5.4) 

For more details about these concepts, see Baccelh and Bremaud [T] or Man- 
nersalo, Norros, and Salminen [Hj. From relationship (j5.3|) we obtain now the 
following property of the joint distribution of Tq and T^. 

Proposition 5.3. For z = or 1 and f > 0, ly > 0, 

P(ri > V, -To >w,Zo = z)=\ P%iTi >t, Zo = z) dt. (5.5) 

J v+w 

Proof: To establish formula ()5.5p we only need to take A = {Ti > v, —Tq > 
w, Zq = z} in fj5.3p and use that for all t G [0, Ti), 

Zqo 9^ = Zq, 

T^odt = T^- t, 

and 

-To oet = t 

P^a.s. (forTo, see dSH). D 

Consider now the on-going busy period of S. There are obvious analogous 
formulas for the idle period. By ()5.5|) . we obtain 

P{-gb>v,dh>w,So>0) = P{T^>v,-To>w,Zo = l) 

= A / P^(Ti >u, Zo = l) du. (5.6) 

Jv+w 
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Consequently, letting 

F;;{v, w) := P{-gb >v,db>w\So>0) 

it is seen that 

F,%v,w) = F,'{v + w,0). (5.7) 

Two-dimensional positive random variables with the property ()5.7|) form a 
special class of random variables studied, e.g., in Salminen and Vallois [T^ . 
The alternative definition of this class, denoted by /C, is as follows. 

Definition 5.4. A two-dimensional positive random variable (X, Y) belongs 

to the class fC if (X, Y) has the same distribution as (f/V, (1 — U)V), where 
the random variable U has the uniform distribution on (0, 1), and V is an 
arbitrary positive random variable independent of U . 

Remark 5.5. In reliability theory the functions with the property (j5.7p are 
called Schur-constant functions (see, e.g., Bassan and Spizzichino [2j). 

The next proposition describes the structure of the joint density of an 
element in /C (for the proof, see Salminen and Vallois 0^). 

Proposition 5.6. Let (X, Y) E IC be such that there exists the density 
Fyiy) =: p{v) . Then 

P(X G dx, Y G dy) = f{x + y) dx dy, 

where 

f{v) = v^-^pi^v), V > 0. 

We have also the following characterization of elements in the class /C in 
terms of the Laplace transforms. 

Proposition 5.7. Let X and Y be positive random variables. Then {X, Y) G 
/C if and only if there exists a positive random variable V such that for all 

E (e-"^-'3^) = ^- r E (e-^^) dj. (5.8) 

o^ ~ P Jp 

In particular, V = X + Y . 
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Proof: =^) Let (X, Y) G /C. Then there exists a positive V such that 

{X,Y)^{UV,{1-U)V), 
where U ~ f/(0, 1) is independent of V. Taking a > P we compute 

= E fe-i»-mv-f3V\ 

E (e-("-/3)^''-/3-) Fvidv) 



1 /"^ e"^^ - e-°^ ^ , , , 



a — P Jq V 



a 



i^/;E(e--).,, 



as claimed in (j5.8j) . 

<^) Now suppose that we have a pair (X, Y) of random variables such that 
()5.8p holds. Letting a ^ P in ()5.8|) and using the continuity of the Laplace 
transforms, we get that X+Y has the same law as V. Further, let U ~ f/(0, 1) 
be independent of V. Then {UV, (1 — U)V) G /C and by the sufficiency part 
of the proof, for a ^ P, 

E ^^-auv-pii-u)v^^ ^ _±^ /"e (e-*^) dt. 

(^ ~ p J p 

Hence, by uniqueness of the Laplace transforms 

(x,r) = {uv,{i~u)v) G/c. 

This completes the proof. D 

Remark 5.8. If {X^Y) G K, and the Laplace transform of X is known, 
we can easily compute the joint Laplace transform. Indeed, denote /(a) : = 
Jp" E (e""^^) d'y. Then setting /5 = in formula ()5.8|) . we get 

E (e--^) = i/(a) 
a 
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and hence, 






a — f3 

We conclude now the proof of Theorem 15. II By ()5.5j) . the two-dimensional 
random variable {—Qbi db) G /C and the marginal Laplace transform of db is 
given by ()3.H) . Hence we can use Remark [5.81 to get the joint distribution of 
—Qf, and db- From Theorem 13. II we have 

* 

E(e-°'=''' I 5*0 > 0) 



r/(a//x)' 
Therefore, the joint Laplace transform of —gb and db is 

a* f a (3 



^^^^g,-m, I ^0 > 0) 



a - /? \Ti{a/fi) r/(/?//i) 



The corresponding formula for {—gi,di) is obtained similarly. The proof of 
Theorem 15. II is now complete. 

6 Example: reflected Brownian motion with 
negative drift 

Let Z = {Zt : t G M} be a reflected Brownian motion with drift — c < in 
stationary state living on J = [0, cx)). Let 

St := supjLj - Ls- {t- s)} 

s<t 

be the local time storage process (as introduced in Section IS} , associated 
with Z. In this case we have 

m{dx) = 2e"^''^'dx, 

and M := m{I} = 1/c. The storage process S is well-defined if and only if 
< c < 1. For a reflected Brownian motion with drift — c the Green function 
at (0, 0) is 

V 2a -t- c^ — c 
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and the function ip (cf. (|2.4|) ) takes the form 

4>{a) = a — v2a + cP + c, a > 

Hence, 

a* = 2(l-c) 

and 

r]{a) = a + l-c+ ^/2aTJl^^c}^ . 

Using ()5.1|) . we compute 
E (e"^''-^'^'' I 5o > 0) 

2(1 - c) / a (5 



a -13 ya + l-c+ v/2a + (l-c)2 /3 + 1 - c + ^2/5 + (1 - c)^ 
1 / 4(1 -c) 4(1 -c) 



a-/3 \v^2/3+(l-c)2 + l + c v^2a + (l-c)2 + l + c^ 
8(1 -c) 



v/2a + (1 - c)2 + ^2/3 + (1 - c)2 

1 



(v/2a + (1 - c)2 + 1 + c)(v/2/? + (1 - c)2 + 1 + c) 
=:F(a,/3;l-c). (6.1^ 

For the on-going idle period at zero, using ()5.2|1 . we obtain after some can- 
cellations that 

E (^gag.-/3d,: I 5^ > 0) =F(a,/?;c). 

This formula can also be found in [9^ but (j6.1|) is a new result. 

To find the joint density of {—gb,db), given that 5*0 > 0, first find the 
density of V^ := —gt + db, given that Sq > 0. Setting /3 = a in the right-hand 
side of ()6.H1 . we get 

E (e-"^'^''-^'') I So > 0) = _^ |^~''^ = . (6.2) 

^ ^ v/2a + (1 - c)2(v/2« + (1 - c)2 + 1 + c)2 

Taking the inverse Laplace transform of ()6.2|) (cf. Erdlyi jH], p. 234) we 
obtain the density of the length of the busy period db — gb, given that 5*0 > 0, 
as 



p,{v) = 2(1 - c) J|^e-^^ - 4(1 - c2)i;e2-$ f-^^\ , v > 0, 
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where $(f ) is the standard normal distribution function. Note that the den- 
sity of the length of the idle period di—gi, given that ^o = 0, ispi_c(w), v > 0. 
Using Proposition Eini we get the joint density of the starting and the ending 
points of the busy period as 

P{—gb G dx, db G dy \ So > 0) 



■K{x + y) 
4(1 - c2)e2'=("'+^)$ /'_i±f ^/^^"j dxdy. 
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